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IDCS J1426+3508: DISCOVERY OF A MASSIVE, IR-SELECTED GALAXY CLUSTER AT Z = 1.75 



(N 

o 

(N 



S. A. Stanford, 1 ' 2 M. Brodwin, 3,4 Anthony H. Gonzalez, 5 Greg Zeimann, 1 Daniel Stern, 6 Arjun Dey, 7 

P. R. ElSENHARDT, 6 GREGORY F. SNYDER, 4 AND C. MANCONE 5 

ApJ, in press 

ABSTRACT 

We report the discovery of an IR-selected massive galaxy cluster in the IRAC Distant Cluster Survey 
(IDCS). We present new data from the Hubble Space Telescope and the W. M. Keck Observatory that 
spectroscopically confirm IDCS J1426+3508 at z = 1.75. Moreover, the cluster is detected in archival 
Chandra data as an extended X-ray source, comprising 54 counts after the removal of point sources. 
We calculate an X-ray luminosity of Lo.5-2/ceV = (5-5 ± 1-2) x 10 44 ergs s _1 within r = 60 arcsec (~ 1 
Mpc diameter), which implies M 20 o,i x = (5.6 ± 1.6) x 10 14 M©. IDCS J1426+3508 appears to be an 
exceptionally massive cluster for its redshift. 

Subject headings: galaxies: clusters: individual — galaxies: distances and redshifts — galaxies: evo- 
lution 
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1. INTRODUCTION 

Galaxy clusters present the opportunity for address- 
ing two main issues in astrophysics. Cosmological pa- 
rameters may be constrained with knowledge of the 
abundance of clusters if the selection function, red- 
shifts, and masses are ad equately known (Haima n et al.l 
1200 It Holder et all 1200 lft . Even individual clusters may 
be useful if they are sufficiently massive and at suffi- 
ciently high redshift because they trace the extreme tail 
of the cosmological de nsity field ()Matarrese et al.l 120001 : 
iMortonson et al.ll2011ft . 

Impressive progress has been made over the past 
decade in finding and characterizing galaxy clusters at 
1 < z < 1.5. At z > 1.5 the number of massive, high red- 
shift clusters that have been identified and confirmed is 
still very limited. The accounting depends on the defini- 
tions both of what constitutes a massive cluster and what 
constitutes confirmation. Here we assume a high-redshift 
cluster massive enough to be useful for constraining cos- 
mological parameters has an M200 mass (the mass within 
the region where the cluster overdensity is 200 times the 
critical density) of at least 1 x 10 14 M©. Such objects 
are the progenitors of present-day clusters with masses of 
i~ 5 x 10 14 M©. We suggest that confirmation of a clus- 
ter candidate requires at least 5 spectroscopic redshift 
members within 2 Mpc. 

At z > 1.5 several systems have been published in 
the literature. The estimat ed masses of th e clust er at 
z = 1.62 identified by both iPapovich et al.l (|2010f ) and 
Tana ka et al.l (12010ft . the cluster at z = 1.75 reported by 
Henry et al.l (J2010), and the cluster at z = 2.07 pub- 
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lished by iGobat et al.1 ( 20~TTI) ap pear to be less than 
1 x 10 14 M©. iSantos et all (|2011l ) and iFassbender et al.l 
(|2011ft have published spectroscopically confirmed X-ray 
selected clusters at z = 1.58 and z = 1.56, respec- 
tively, both of which appear to have cluster masses of a 
fewxlO 14 M©, although in the former case only 3 mem- 
ber galaxies have spectroscopic redshifts. Thus there are 
only 1-2 galaxy clusters at z > 1.5 that are massive, con- 
firmed, and published. A number of prot o-clusters have 
been identified and confirmed at z > 2 (TPentericci et al.1 
120001 IVenemans et al.l 120071: ICanak et all 12011ft . While 
very interesting, these systems appear to be in the very 
early stages of cluster formation-their masses remain dif- 
ficult to estimate and their natures difficult to interpret. 

The other main astrophysical use for galaxy clusters 
is to help us understand galaxy evolution. In particu- 
lar, they contain the majority of the massive early-type 
galaxies in the universe, so these environments offer ex- 
cellent places in which to study such galaxies. To trace 
the evolution of massive early-type galaxies over their 
full lifetime, we must identify and study the precursor 
cluster population over a large redshift range. This kind 
of archaeology requires the evolutionary precursors to be 
identified in large, statistically useful samples, which are 
sensitive down to the group scale at relatively high red- 
shift, since the massive clusters at z < 1 are built from 
groups and low-mass clusters at z > 1. 

Identifying and then characterizing cluster samples ad- 
equate to these two tasks has been challenging. Op- 
tical methods of finding clusters tend to succeed up 
to redshifts only slightly beyond unity. Neither X-ray 
nor Sunyaev-Zel'dovich (SZ) cluster surveys currently 
have the sensitivity to reach cluster masses down to 
1 x 10 14 M© at z ~ 1.5 and above over appreciable areas. 

The IRAC Shallo w Cluster Survey (ISCS, 
lEisenhardt et al.1 12008ft originally was designed to 
create a stellar mass-selected selected sample of galaxy 
clusters spanning < z < 2. The IS CS is drawn from 
the Spitzer/IRAC Shallow Survey ()Eisenhardt et al.l 
2004), which imaged most of the Bootes field in the 
NOAO Deep Wide-Field Survey (|.Tannuzi fc Devi fl999t) . 
Clusters were identified by searching for 3-dimensional 
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spatial overdensities in a 4.5/im s elected galaxy sampl e 
with robust photometric redshifts (|Brodwin et al.ll2006h ■ 
The selection is independent of the presence of a red 
sequence. There are 335 clusters and groups in the ISCS 
sample, identified over 7.25 deg 2 within the Bootes field, 
and 1/3 of the groups/clusters are at z > 1. We have 
spectroscopic ally confirmed over 20 clusters spanning 
1 < z < 1.5. ([Stanford et al.l 120051: IBrodwin et aL|[200l 
Elston et aLll2006l : lEisenhardt et alj|2008t iBrodwin et al.l 
20111 Brodwin et al. in prep, Zeimann et al. in prep). 

We have begun extending the work of the ISCS to more 
effectively target the higher redshift range by making 
use of the deeper IRAC imaging obt ained by the Spitze r 
Deep, Wide Field Survey (SDWFS. lAshbv et al.l[2009h . 
The IRAC exposure time of SDWFS is 4x that of the 
IRAC data used in the ISCS, which has allowed us to bet- 
ter identify cluster candidates at z > 1.5 in the Bootes 
field as part of the new IRAC Distant Cluster Survey 
(IDCS). Here we report the first spectroscopic confirma- 
tion of one of these candidates, IDCS J1426+3508 at 
z = 1.75, which we believe to be a massive cluster. Com- 
panion papers in this volume report on the mass of this 
cluster (Brodwin et al. 2012) and on a gravitational arc 
(Gonzalez et al. 2012). We present the optical and in- 
frared imaging in |J21 the spectroscopic observations and 
resulting redshifts in and [JU respectively, and the X- 
ray observations in JjS We use Vega m agnitudes and a 
WMA P7+BAO+ff ACDM cosmology (jKomatsu et al.1 
[20T1 : fl M = 0.272, fl A = 0.728, and H = 70 .4 km s" 1 
Mpc^ 1 . 

2. OPTICAL AND NEAR-IR IMAGING 

The cluster candidate was originally identified using 
the SDWFS data matched with the NDWFS optical 
da ta following the same p rocedures described in detail 
in lEisenhardt et al.1 (j2008j ) . A significant overdensity in 
the 3-dimensional space of (RA, Dec, and photometric 
redshift) was selected for further study. A color image 
made from the NDWFS optical + IRAC data is shown 
in the left panel of Figure [T] where a tight red group of 
galaxies is visible. The photometric redshift estimate of 
the cluster is z ~ 1.8. 

Deeper optical and NIR follow-up imaging was ob- 
tained with the Hubble Space Telescope using ACS and 
WFC3. F814W exposures were obtained with ACS in 
one pointing for 8 x 564 s and reduced using standard 
procedures. The WFC3 data were obtained with the 
F160W filter in two slightly overlapping pointings, each 
comprising 700 s of integration time in dithered expo- 
sures. The WFC3 data were reduced using standard pro- 
cedures with the MultiDrizzle software. A pseudo-color 
image constructed from the registered ACS and WFC3 
imaging is shown in the right panel of Figure [TJ 

Photometry w as measured from the H ST imaging us- 
ing SExtractor (jBertin fc Arnouti Il996) in dual image 
mode with sources detected in the WFC3 image. Col- 
ors were measured in 0.8 arcsec diameter apertures, and 
MACAUTO is used as a "total" magnitude. Our pho- 
tometric uncertainties are dominated by sky shot noise. 
We estimated the uncertainties from the distribution of 
sky background measurements in 5000 randomly placed 
0.8" apertures. These sky measurements are roughly nor- 
mally distributed; we estimate a s k y by fitting a Gaussian 
to the left-half of the distribution (which is uncontami- 



nated by light from objects). We verified this procedure 
by confirming that it produces the correct scaling in pho- 
tometric scatter of sources detected in the sets of dither 
images before making the final stack. 

Morphologies of the galaxies i n the WFC3 image w ere 
determined by using Galapagos ( Hausslcr et al.ll20TlT ) to 
run GALFIT (|Peng et al.l f2010h . Galapa gos measures 
the sky around every galaxy and uses the basic isopho- 
tal parameters measured by source extractor to generate 
a first guess set of Sersic parameters for each galaxy. 
It then uses these values to fit a single Sersig profile to 
every galaxy with GALFIT, simultaneously fitting close 
neighbors. We use the Sersig index n to classify galaxy 
morphologies as being either early-type (n > 2.5) or late- 
type (n < 2.5). 

3. SPECTROSCOPY 

3.1. Keck Optical Spectra 

Spectra wi th the Low-Res olution Imaging Spectro- 
graph (LRIS; lOke et afl 119951 ) on Keck I were acquired 
on UT 2011 April 28 and 29 using 1.1" x 10" slitlets, the 
G400/8500 grating on the red side, the D680 dichroic, 
and the 400/3400 grism on the blue side. On the 
first night when the seeing was 0'.'6 and conditions were 
mostly clear, data were obtained in one mask in six 1200 s 
exposures. On the second night when the seeing var- 
ied in the range 0.6 — 0.8 arcsec with light cirrus, data 
were obtained in a second mask for ten 1200 s expo- 
sures. These data were split into slitlets which were sep- 
arately reduced following standard procedures. The rela- 
tive spectral response was calibrated via longslit (with a 
wid th of 1.0 arcsec) observati ons of Wolf 1346 and Feige 
34 (|Massev fc GronwaHl[l990( ) . 

3.2. WFC3 Grism Spectra 

The WFC3 data were obtained on 2010 November 6. 
A total exposure time of 11247 s was used for the G102 
observation and 2011 s for G141; a single pointing was 
used in both cases. We targeted the cluster candidate 
with both IR prisms so as to provide continuous wave- 
length coverage of nearly 1 /im, allowing for conclusive 
identifications of spectral features in the redshift range 
of interest. The FWHM of an unresolved emission line is 
^2 pixels, which corresponds to a spectral resolution of 
49 A for G102 and 93 A for GUI. The low grism resolu- 
tion blends some of the more common emission lines such 
as Ha+N[II] and [OIII]5007+[OIII]4959. This typically 
results in a redshift uncertainty of a z ~ 0.01. 

In slitless spectroscopy, a direct image is a necessary 
companion to the grism image in order to zero-point 
the wavelength scale and to properly extract spectra. 
We chose broadband filters that closely matched the 
grism spectral coverage: F105W for G102 and F140W 
for GUI. The object positions and sizes measured from 
the direct images are used to establish the location, wave- 
length zero point, and spectral extraction widths of the 
objects in grism images. The data were reduced using 
aXe (Kummel et al. 2008). The calibration filefQ used 
were the best available at the time of the reductions. 
The steps used to extract spectra are very similar to that 
found in WFC3 Grism Cookbook and a more detailed 

1 http: / / www.stsci.edu/hst / wfc3 / analysis / grism_obs / calibrations / 

2 http: / / www.stsci.edu/hst /wfc3/analysis / grism_obscookbook.html 
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Fig. 1. — (left) Color image covering 3x3 arcmin using imaging from the NDWFS Bw and /, and IRAC 4.5 /an data centered on IDCS 
J1426+3508. The Q marks the quasar in the cluster, (right) Pseudo-color HST image made from the ACS/F814W and WFC3/F160W 
images with the green contours illustrating the X-ray emission. The dashed red circle is centered on the quasar in the cluster, and the blue 
dashed circle is centered on a non-member radio-loud AGN. The radii of these two dashed circles is 5 arcsec, the same size as was used to 
mask these point sources in the X-ray analysis. In both panels the yellow boxes are spec-z confirmed members, and a 30" (260 kpc) scale 
bar is given. 



description will be presented in Zeimann et al. (in prepa- 
ration) . 

4. REDSHIFT MEASUREMENTS 

The reduced optical spectra were visually inspected to 
determine redshifts. Despite the long integrations on 
a 10 m telescope using a spectrograph with new red- 
sensitive CCDs, sufficiently good LRIS spectra were ob- 
tained on only two phot-z selected objects which were 
found to have features such as D4000, B2640 and the 
MgIIA2800 absorption line that are characteristic of older 
stellar populations. These features indicate z = 1.75 for 
both objects, which were confirmed by the WFC3 IR 
grism spectroscopy. 

The NIR grism spectra were first visually in- 
spected. Emission lines were identified as be- 
ing blended Ha+[NII], [OIII]A5007+[OIII] A4959, 
H/3, or [OII]A3727. The spectra were also cross- 
correlated with spectral templates (taken from 
www.sdss.org/ dr7/algorithms/spectemplates/index.html) 
to automatically determine redshifts where feasible 
(Zeimann et al., in prep). Contamination from overlap- 
ping spectra was estimated using a Gaussian emission 
model, scaled by the measured broadband magnitudes 
from the direct images (F105W or F140W; see the 
WFC3 Grism Cookbook for more details). This process 
is handled in the standard reduction package of aXe. 
For the ETG spectra, we experimented with a range 
of SDSS galaxy templates and found that the LRG 
template was the best match to our grism data. 

The spectra of cluster members are shown in Figure [5] 
and Figure [3] The latter shows in the bottom panel 
a simulated early-type galaxy spectrum, as observed by 
WFC3 with our observational parameters. aXeSim was 



used with an SDSS LRG template spectrum^, redshifted 
to z = 1.75, to create a mock grism image. A 1-D extrac- 
tion was performed with the same reduction procedure 
as was used with the actual WFC3 grism observations. 
The redshifts determined from both the LRIS and WFC3 
spectroscopy on the identified cluster members are sum- 
marized in Table [TJ along with magnitudes and colors 
obtained from the HST images. The overall redshifts ob- 
tained in the vicinity of the cluster candidate with the 
LRIS and WFC3 grism data are presented in Figure |H 

5. X-RAY OBSERVATIONS 

The Bootes field has been surveyed previously 
with ACIS-I onboard the Chandra X-Ray Observa- 
tory (iMurrav et al.ll2005l : iKenter et ai1l2005t iBrand et al] 
l2006h . At the position of IDCS J1426+3508, exposures 
totalling 9.5 ks are available from the Chandra archive. 
These data are split between an observation of 4.8 ks on 
UT 2006 July 30 (ObsID 3621) and an observation of 4.7 
ks on UT 2006 August 21 (ObsID 7381). We processed 
the data following standard procedures using the Chan- 
dra Interactive Analysis of Observations (CIAO; V4.2) 
software. We initially identified good-time intervals for 
the exposures, yielding a total effective exposure time of 
8.3 ks for IDCS J1426+3508. 

The cluster is clearly detected as an extended source 
in both individual exposures, as well as in the stacked 
exposure. The cluster is approximately 6.5 arcmin off- 
axis in both exposures, for which the Chandra point- 
source 50% encircled energy radius is 2 arcsec at 1.5 
keV. This complicates the X-ray analysis, as an opti- 
cally bright quasar confirmed to be in the cluster (see 

3 http:/ /www. sdss.org/dr7/algorithms/spcctcmplatcs/indcx. html 
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Fig. 2. — WFC3 spectra of the five cluster members that exhibit emission lines are plotted above. The solid black histograms are the 
spectra from the G102 and G141 grisms. The dot-dashed cyan line is the estimate of contamination from overlapping spectra which is 
subtracted off in the final stage of reduction. The solid red line is the 1-cr flux error. The vertical green lines which are labeled are the 
detected or expected emission from the [OIIJA3727, H/3, and [OI11]A5007 lines at the nominal cluster redshift. The bright, power-law 
spectrum second from the top is a QSO, previously identified in AGES optical spectroscopy (Kochanek et al., in preparation); in this panel 
a QSO template (SDSS; Vanden Berk et al. 2001) is shown by the magenta line. 



TABLE 1 

Spectroscopic Cluster Members 



ID 


R.A. (J2000) 


Dccl. (J2000) 


2 


Az 


Instrument 


F160W 


F814W-F160W 


J142632.9+350823 a 


14:26:32.95 


35 


08 


23.6 


1.75 


0.01 


WFC3/LRIS b 


19.25 


4.47 


J142632. 5+350822 


14:26:32.55 


35 


08 


22.5 


1.75 


0.01 


WFC3/LR1S C 


20.41 


4.16 


J142632.4+350830 d 


14:26:32.40 


35 


08 


30.8 


1.746 6 


0.01 


WFC3 


18.59 


1.08 


J142632. 8+350844 


14:26:32.85 


35 


08 


44.4 


1.74 


0.01 


WFC3 


22.81 


1.20 


J142634.4+350825 


14:26:34.43 


35 


08 


25.1 


1.75 


0.01 


WFC3 


22.63 


1.14 


J142634.4+350822 


14:26:34.47 


35 


08 


22.4 


1.75 


0.01 


WFC3 


22.27 


0.90 


J142628. 1+350829 


14:26:28.15 


35 


08 


29.7 


1.75 


0.01 


WFC3 


23.00 


1.43 



a Brightest Cluster Galaxy 

LRIS spectrum shows a break at 2640 A and a very red continuum consistent with z — 1.75 
c LRIS spectrum shows a red continuum and a MgIIA2S00 absorption feature consistent with z — 1.75 
d QSO 

e AGES redshift; WFC3 grism redshift is 1.77 
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Fig. 3. — WFC3 spectra of the two cluster members with early-type spectra in the top two panels; the bottom panel shows a simulated 
spectrum, described in the text, for reference. The solid black histograms are the spectra from the G102 and G141 grisms. The dot-dashed 
blue line is the estimation of contamination from overlapping spectra which is subtracted off in the final stage of reduction. The solid red 
line is the 1-ct flux error. The vertical green lines are the expected locations of the [OIIJA3727, H/3, and [OIII]A5007 lines. The vertical 
red lines are the expected locations for the following absorption features: Ca H+K, the G-band, H7, and MgIIA2800. The magenta lines 
represent the SDSS LRG template fitted to the observed spectra. The inset spectrum in the top panel shows the LRIS spectrum (black 
solid line) in the vicinity of the D4000 break, along with the template fit (magenta solid line), solidifying the reality of this feature seen in 
the grism spectrum. 
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Fig. 4. — Redshift histogram resulting from the LRIS and WFC3 
spectroscopic observations. The green bars show all spectroscopic 
redshifts of all qualities. The blue histogram shows only the robust 
redshifts (i.e., quality A or B), and the red histogram shows the 
subset of these which lie within 30" of the brightest cluster galaxy. 
The inset shows a detail of the redshift histogram near the cluster 
redshift; its width is consistent with the uncertainties of the grism 
redshifts. 

Section 5) is only 9 arcsec from the brightest cluster 
galaxy (BCG). In addition, we identify an X-ray point 
source associated with a radio source that is 12 arcsec 
to the SW, at 14:26:32.2, 4-35:08:14.9. This source, con- 
firmed as an emission line galaxy at z = 1.535 in our 
WFC3 grism data, has an integrate d 21 cm flux densit y 
of 95.3 mJy from the FIRST survey (jBecker et al.lll995| ). 
Given the signal to noise ratio and the large off-axis an- 
gle of the available X-ray observations, it is challenging 
to disentangle the extended cluster emission from the 
point source contributions. However, as seen in the right 
panel of Figure [Q IDCS J1426+3508 is clearly associ- 
ated with diffuse X-ray emission that extends beyond 
the point-source contributions from the two AGN. We 
see no evidence in the WFC3 grism spectroscopy, which 
covers the central 2 arcmin of the cluster, for other AGN 
which could contribute to the X-ray emission. We have 
also used the SDWFS IRAC photometry to construct a 
two-color diagram to search for obscured AGN (following 
Stern et al. 2005) which might contribute to the mea- 
sured X-ray flux. In addition to the QSO and the radio 
source already described, there is one more object within 
the X-ray measurement aperture which has IRAC colors 
typical of obscured AGN. 

To extract the X-ray counts due to the cluster, we 
masked the three AGN (the QSO, the radio source, 
and the IRAC AGN) using a conservative 5 arcsec ra- 
dius aperture, corresponding to the 90% encircled en- 
ergy radius at the observed off-axis-angle of the clus- 
ter. We expect that the unmasked flux from the two 
AGN to contribute only one photon to the measured 



flux. We then extracted cluster source counts in the 
0.5 - 7 keV range within a 1 arcmin radius aperture 
centered on the cluster BCG. This aperture approxi- 
mately corresponds to a radius of 500 kpc at the clus- 
ter redshift. Response matrices and effective areas were 
then determined for each detected source. Within the 
measurement aperture, there are 54 ± 10 background- 
corrected counts in the 0.5 — 7 keV range, after masking 
out the two central AGN. We used XSPEC (V12.6.0) to 
fit the background-subtracted X- ray spectrum with the 
MEKAL hot, diffuse gas model (jMewe et all 119851 ) us- 
ing the Wisconsin photo-electric absorption cross-section 
(|Morrison &: McCammonl 119831 ). The temperature was 
fixed at 5 keV and the abundance at 0.3 M©, with a 
Galactic absorption of 1.3 xl0 20 cm~ 2 at the target po- 
sition. We determined a Galactic absorption-corrected 
flux of (3.1 ± 0.7) x 10" 14 ergs cm" 2 s" 1 in the 0.5 - 2.0 
keV range, which translates to an X-ray luminosity of 
(5.5 ± 1.2) x 10 44 ergs s" 1 at z = 1.75. The X-ray flux 
changes by only 7% if the X-ray temperature is varied 
from 4 to 6 keV. 

Using the M 50u -Lx relation of lVikhlinin et al.l (|2009l> . 
we estimate from the luminosity that Msoo,l x = (3.5 ± 
1.0) x 10 14 Mq. We caution that the use of this scal- 
ing relation requires a significant extrapolation in red- 
shift. This and other systematic uncertainties (such 
as removal of X-ray point sources) are expected to 
dominate over statistical errors. To estimate a total 
clus ter mass, we next assume an NFW density pro- 
file (INavarro et al.lll997j) and the mass-concentration re- 
lation of (jDuffv et al.l (20081) . The resulting M 200 = 
(5.6 ± 1.6) x 10 14 Mq, where the uncertainty is deter- 
mined only from the statistical measurement error. 

6. GALAXY POPULATIONS 

Figure O presents the photometric and morphological 
information for all the objects in the 7.7 arcmin 2 area 
centered on the cluster where the ACS and WFC3 imag- 
ing overlap. While it is possible to see a red sequence of 
early-type galaxies in Figure[5l the spread and location is 
different from that of the red sequence in massive clusters 
at 1 < z < 1.5. To isolate a sample of probable cluster 
members in the color-magnitude diagram we carry out 
the following multi-step procedure. First, we determined 
the relative offset A between the object colors and the 
expected color-magnitude relation at t his redshift (based 
on a m odel for Coma with Zf = 3; see [Eiscnhard t et al.l 
([2001 for details). The range -0.5 < A < 1.0 (corre- 
sponding to colors 2.3 < F814W - F160W < 3.8) is the 
color cut used to select member galaxies on the red se- 
quence. We initially define the cluster red sequence as 
galaxies brighter than H * (z) + 1 .5 that are chosen by the 
above color cut in A, where H*(z) is passively evolved 
from the luminosity of Coma early- type galaxies. Most 
of the objects in this color-selected red sequence do not 
have spectroscopic or photometric redshift information, 
and so the initial red sequence sample may suffer from in- 
terloper contamination. For the purpose of studying the 
cluster red sequence, we choose to restrict consideration 
to morphologically selected early-type galaxies, keeping 
objects with n > 2.5. We discard objects with colors 
more than two absolute deviations from the central red 
sequence color, where the deviation is the median of the 
absolute value of the A of the potential red sequence ob- 
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jects. Objects removed in the latter step are shown with 
a gray dot in Figure [6] The remaining objects in the 
red sequence sample are represented by the red points in 
Figure El 

To measure the color and scatter of this red sequence 
sample down to F160W = 21.8 6, we use the biw eight 
estimates of location and scale ([Beers et al.l lT990). We 
calculate the intrinsic scatter Oi n t in the red sequence 
sample by subtracting in quadrature the median color 
error from the biweight scale estimate. The red hori- 
zontal lines in Figure |6] arc offset a distance <7j n t above 
and below the median color of the red sequence. We 
estimate uncertainties by performing these calculations 
on 1000 bootstrap resamplings, from which we measure 
the scale of the resulting color and (Ji n t distributions us- 
ing their median absolute deviation. The intrinsic scat- 
ter in the observed colors of the red sequence galaxies is 
0.16 ±0.06 mag; transformed to the rest frame U — V the 
intrinsic scatter is 0.10 ± 0.04. The median color of the 
morphologically-selected red sequence sample is F814W 
- F160W = 3.96 ±0.07, which is 0.29 ±0.05 redder than 
the expected color, derived fr om a simple passiye evo- 
lution model with a Zf — 3 (|Bruzual fc Charlotl 120031 
2007 version). A model with an earlier formation epoch, 
Zf = 6, predicts a color which is closer to the observed 
F814W — F160W colors of the morphologically-selected 
red sequence. 

The cluster includes a fairly large brightest cluster 
galaxy. Profile fitting using GALFIT determined a Sersig 
index of 5.4 ± 0.1, r e = 18.0 ± 2.5 kpc, and a to- 
tal F160W magnitude = 18.5 ± 0.05 which is approxi- 
mately 2.0 magnitudes brighter than M* for a passively- 
evolv ing galaxy formed at Zf = 3 (|Bruzual fc Charlotl 
120031 ). From the F160W total magnitude, we determine 
a rest frame My = —24.7 assuming the Zf = 3 passive 
evolution model to transform to rest V-band. The size 
and luminosity of this BCG is similar toaBCGatz~0, 
which is remarkable since it is at z = 1.75. 

7. DISCUSSION 

IDCS J1426±3508 is a newly discovered galaxy clus- 
ter at z — 1.75, which places it among the few such 
systems currently known at 1.5 < z < 2, when mas- 
sive clusters may be first forming. Using a combination 
of optical multi-object Keck spectroscopy and infrared 
HST/WFC3 grism spectroscopy, we have confirmed 7 
cluster members in IDCS J1426+3508 within a radius 
of 2 Mpc, all but one of which are within a radius of 
250 kpc. The extended X-ray emission described in Sec- 
tion 4 indicates that this cluster is gravitationally bound 
and already has a mass greater than 10 14 M©. 

The properties of the galaxies in the cluster indicate 
that the cluster is far from settled in terms of star forma- 
tion. Most of the spectroscopic member galaxies are very 
blue (c.f. Table [IJ and show emission lines in their spec- 
tra indicative of on-going star formation or AGN activity. 
The red sequence itself has a larger amount of scatter, 
0.16 mag, compared to 0.08 in a similar observed color 
in XMMU J2235— 2557, a slightly mor e massive cluster 
at so mewhat lower redshift, z — 1.39 (jStrazzullo et al.l 
12010ft . The intrinsic scatter in IDCS J1426±3508 is rea- 
sonable for a population of galaxies whose stars formed at 
Zf ~ 3.5 — 6. The median color of these galaxies matches 
the prediction for simple passive evolution if the stars in 



the red sequence galaxie s were formed at 5 < Zf < 6 
(|Bruzual fc Charlotl 120031 2007 version). This is in con- 
trast to the cluster at z = 1.62 (jPapovich et al.l 120101 : 
iTanaka et al.l [2010) in which the red sequence has a 
younger Zf. The red sequences may have different me- 
dian colors because of the difference in the cluster masses, 
in that IDCS J1426±3508 is more massive and so may 
have formed earlier. Reddening by dust could cause the 
redder colors in IDCS J1426±3508, though there is no 
reason to suspect unusual amounts of dust in this cluster. 
Indeed, most of the spectroscopically confirmed members 
are fairly blue. 

8. CONCLUSIONS 

We have presented optical and NIR spectroscopy of 
galaxies in IDCS J1426±3508, an IRAC-selected galaxy 
cluster candidate, which indicate this is a bona fide clus- 
ter of galaxies at z = 1.75. The available Chandra data 
show a faint but clearly extended X-ray source at this 
location. Along with the centrally concentrated and reg- 
ular distribution of red galaxies, the X-ray detection ar- 
gues that IDCS J1426±3508 is gravitationally bound and 
in a relatively relaxed state. 

The discovery of a cluster at z = 1.75 with a soft X- 
ray flux greater than the ~ 10 -14 e rgs cm" 2 s" 1 limit o f 
the upcoming eROSITA mission (IPredehl et al.ll20Tol) . 
scheduled for launch in 2014, promises exciting results 
for that mission. However, we note that at the 28 arc- 
sec average resolution expected for the eROSITA slew 
survey, and even at the 15 arcsec on-axis resolution of 
eROSITA, it would have been very difficult to distin- 
guish the AGN from the cluster emission. Without the 
AGN masked out in our Chandra data, the derived X- 
ray flux for IPCS J1426±3508 would be - 50% higher. 
As shown by e.g. iGalametz et al.l (|2009ft . the fraction of 
galaxy clusters that host luminous AGN increases rapidly 
with redshift. Such AGN will make it challenging to draw 
firm cosmological evolutionary results from the eROSITA 
cluster sample without higher resolution X-ray follow-up. 

With the higher resolution Chandra data currently 
available from the archive, the point-source corrected X- 
ray luminosity L0.5-2.0keV — (5.5 ± 1.2) x 10 44 ergs s _1 , 
which implies M2oo,l x ~5x 10 14 M Q . This is a surpris- 
ingly large cluster mass for this redshift and survey area, 
as will be addressed in forthcoming papers. 

This work is based in part on observations obtained 
with the Chandra X-ray Observatory (CXO), under con- 
tract SV4-74018, A31 with the Smithsonian Astrophys- 
ical Observatory which operates the CXO for NASA. 
Support for this research was provided by NASA grant 
G09-0150A. AHG acknowledges support from the Na- 
tional Science Foundation through grant AST-0708490. 
This work is based in part on observations made with 
the Spitzer Space Telescope, which is operated by the 
Jet Propulsion Laboratory, California Institute of Tech- 
nology under a contract with NASA. Support for this 
work was provided by NASA through an award issued 
by JPL/Caltech. Support for HST programs 11663 and 
12203 were provided by NASA through a grant from 
the Space Telescope Science Institute, which is oper- 
ated by the Association of Universities for Research in 
Astronomy, Inc., under NASA contract NAS 5-26555. 
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19 20 21 22 

F160W 

Fig. 5. — Color-magnitude diagram of IDCS J1426+3508, made from ACS and WFC3 imaging. The red and blue colored points denote 
morphological classifications, based on the Sersig index. The spectroscopically confirmed members are marked by the larger stars; only 
two of the members are red enough to appear in this plot. The black dotted line repr esents the expected col or of a passively evolving red 
sequence of galaxies formed at zj = 3, with the slope based on observed Coma colors (Eiscnhardt ct al. 2007). The diagonal black dashed 
line is the fit to the red sequence galaxies. The diagonal dashed blue line represents the 5 a limit on the colors. 




F160W 

Fig. 6. — The left panel of the plot is similar to the color-magnitude diagram shown in Figure[5] The F814W — F160W colors have been 
zero pointed to the color predicted by a passively evolving red sequence of galaxies formed at Zf = 3, with the slope based on observed 
Coma colors (Eiscnhardt ct al. 2007). The red filled points are the objects selected by the magnitude and color cuts described in the text 
which have a Sersig index n > 2.5. The grey-filled circles were removed from the red sequence selection. The solid red horizontal lines 
show the one (Ti„t intrinsic scatter of the red sequence galaxies, and the black dashed line is the fit to the colors of these galaxies. The two 
spectroscopically confirmed members that lie within the limited color range of this CMD are marked by the larger stars. The right side of 
the plot is a histogram stacked in the colors. 



IDCS J1426+3508 at z = 1.75 
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